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Abstract 



We study higher-order variational derivatives of a generic second-order Lagran- 
gian C = C(x, (ft, deft, d 2 (ft) and in this context we discuss the Jacobi equation ensuing 
from the second variation of the action. We exhibit the different integrations by 
parts which may be performed to obtain the Jacobi equation and we show that 
there is a particular integration by parts which is invariant. We introduce two new 
Lagrangians, C\ and £2, associated to the first and second-order deformations of 
the original Lagrangian Co respectively; they are in fact the first elements of a whole 
hierarchy of Lagrangians derived from Cq. In terms of these Lagrangians we are able 
to establish simple relations between the variational derivatives of different orders of 
a given Lagrangian. We then show that the Jacobi equations of £$ may be obtained 
as variational equations, so that the Euler-Lagrange and the Jacobi equations are 
obtained from a single variational principle based on the first-order variation C\ of 
the Lagrangian. We can furthermore introduce an associated energy-momentum 
tensor Tt^ u which turns out to be a conserved quantity if Cq is independent of 
space-time variables. 

Introduction 

As is well known, the second variation of an action functional governs the behaviour of 
the action itself in the neighbourhood of critical sections. In particular, the Hessian of the 
Lagrangian defines a quadratic form whose sign properties allow to distinguish between 
minima, maxima and degenerate critical sections (see, e.g., [1]). It is also well known that 
in the case of geodesies in a Riemannian manifold those fields which govern the transition 
from geodesies to geodesies, (i.e., those vectorfields which make the second variation to 
vanish identically modulo boundary terms) are Jacobi fields [2] and they are solutions of 
a second-order differential equation, namely the Jacobi equation (of geodesies). 

The notion of Jacobi equation as an outcome of the second variation is in fact fairly 
more general than this and general formulae for the second variation and generalised 
Jacobi equations along critical sections have been already considered in the Calculus of 
Variations from a "structural viewpoint" (see, e.g., the review of results contained in [3]). 
To our knowledge, however, in most of the current literature on the subject the second 
variation of functionals has been considered in a direct way and without resorting to 
general expressions, while integration by parts to reduce it to more suitable forms have 
been usually performed by ad hoc procedures, in spite of the fact that fairly general such 
formulae can be worked out (see below for a comment on a paper by Taub [4], which is 
appropriate to mention here but is better to discuss later). Because of these facts, we have 
reached the conclusion that the theory of second variations is worthy of being revisited, 
also in view of a number of applications which shall be mentioned later and will form the 
subject of forthcoming investigations. 

In a previous paper [5], working in the general and global framework of first-order 
variational principles in fibered manifolds and their jet-prolongations, the notion of gen- 
eralised Jacobi equations was developed and discussed for first-order Lagrangians. As 
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is well known, the classical Jacobi equation for geodesies of a Riemannian manifold de- 
fines in fact the Riemannian curvature tensor of g; because of this we can say that the 
second variaton S 2 S and the generalised Jacobi equations define the "curvature" of any 
given variational principle. In the generic case, of course, this notion has very little to 
say, in a continuation [6] of that work it was however shown that this general concept of 
"curvature" takes a particularly significant form in the case of generalised harmonic La- 
grangians, giving rise to suitable "curvature tensors" which satisfy suitable "generalised 
Bianchi identities". Applications to second variations of relativistic Lagrangians (i.e., La- 
grangians depending on the full curvature of a Riemannian metric) are being considered 
in [7] and [8]. 

For the sake of completeness, we should mention that a number of recent and less 
recent papers (see, e.g., [9], [10] and references quoted therein) have attacked the problem 
of "generalised Jacobi equations" and "curvature" for arbitrary second-order differential 
equations, from the general viewpoint of dynamical systems on a tangent bundle. These 
interesting contributions have a rather different nature from ours, since they refer to 
fairly general dynamical structures while the results we are going to discuss stem directly 
from the richer structure of Lagrangian systems on generic fiber bundles. Nevertheless, 
it would be interesting to compare the two viewpoints, also to check how much of our 
direct, simpler and straightforward results might be recovered from a suitable application 
of the more complicated structures discussed in [10] (and references quoted therein) in 
a different and complementary framework. We hope to address this problem in future 
investigations. 

In this paper we shall consider some further results related to the generalised Jacobi 
equation. For the sake of simplicity and for mere notational convenience we will estab- 
lish our results only for second-order Lagrangians in field theory, although an analogous 
scheme can be worked out for higher-orders without bringing any really new insight into 
the problem. 

The first part of this paper (Section 1) is devoted to briefly recall the situation for 
first-order theory, as already discussed in [11] in the intrinsic language of differential ge- 
ometry of jet-bundles. In Section 2 we consider the second variation 5 2 S of the action S, 
defined by a Lagrangian depending at most on second-order derivatives of the fields. We 
then show how a number of different integrations by parts allow to recast the Hessian in 
more suitable forms, which contain the Euler-Lagrange equations and define some ordi- 
nary differential equations of the second-order which are the generalised Jacobi equations. 
These are in fact the equations which define along critical curves those vectorfields in the 
configuration space which make 5 2 S to vanish identically (modulo boundary conditions). 

In Section 3 we introduce then a series of relevant relations between the Euler- 
Lagrange equation for a given Lagrangian C , namely 

5 _£± = f?£o _ J_ ( dC \ d 2 ( dC \ 

54> A d<p A dx»\d(j) A J dxudxnyd^^J ' l ' j 

and the "Jacobi equation" for the same Lagrangian £ , namely 
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d 2 C 
d 2 C 

n'\ • , ff" 'I : 0. (0.2) 




dx"dx" \d^ d(f) B I d^ u dcf> B x ' d$%d<t> B Xp 



Starting from the Lagrangian C we shall in fact define new Lagrangians C\ and £ 2 
associated respectively to the first-order deformation of Co and to the second-order de- 
formation of Co, with the remarkable property that the Euler-lagrange equations (0.1) 
and the Jacobi equations (0.2) of the original Lagrangian Co can be rewritten together in 
terms of C\ and C 2 as follows: 



SCq 5C\ 

S(f> A Sr] A 
5C 2 _ 6C1 
Srj A S(j) A 



0, (0.3) 
0. (0.4) 



Therefore, we can consider the above relations as coming from a single variational principle 
based on the new Lagrangian C\ where, however, 0's and 77's are considered as independent 
variables. This result finds important applications to Riemannian Geometry [11]. 

We can thence introduce momenta canonically conjugated to 0's and 77's by means 
of standard prescriptions (see Section 3 below) as well as a canonical energy-momentum 
tensor 7"P„, which is conserved if Co does not depend explicitly on "space-time variables" 
x x . Therefore, the allowed deformations are selected by the first-order condition 

—^ = 0, 0.5 

which replaces the second-order condition given by the Jacobi equation (0.2). 

We are now in position to mention in a greater detail the paper by Taub [4] we already 
quoted above. In that paper, which is explicitly devoted to an application of second- 
order variations to relativistic fluid-dynamics, a general formula for the Jacobi equation 
is explicitly mentioned and it is claimed that a "well-known result" states that this 
equation is in fact the variational derivative of the Hessian C 2 (see equation (0.4) above. 
Unfortunately, no reference is given to the source of this result; we stress, however, that 
the viewpoint we shall discuss in this paper is fairly different and more general. In fact, 
we shall not see Jacobi equations alone as (partial with respect to 77) Lagrange equations 
of the Hessian C 2 (as in [4]), but rather the system formed by Jacobi equations and the 
original Lagrange equations of Co as variational equations of the first-order deformed 
Lagrangian C\. 

Section 4 is finally devoted to some applications and examples. 
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For the sake of simplicity most of this paper will be written in the language of fun- 
damental calculus in R™, although all the results presented here can be expressed and 
derived in the intrinsic language of fiber bundles. 

1 The Lagrangian Characterization of Generalized 
Jacobi Equation in the First-Order Case 

Let us first recall some basic notion from the Calculus of Variations on fibered manifolds. 
Let (B, M, 7r) be a fibered manifold and £ : J 1 ^ A n (T*M) be a fisrt-order Lagrangian 
(density) on B. Here J l n denotes the first-jet prolongation of 7r and A n (T*M) is the 
bundle of n-forms of M, n being the dimension of the base manifold M. Locally 

C = L(x\ <f> A , 4> A x )ds, (1.1) 

where ( J 1 £7; x x , <p A , 4> A \) is any natural chart in J lr K and ds is the local volume of (U ; x x ). 
The action of C is defined by 

S = ( (fa)* C , (1.2) 
Jn 

where Q is any compact domain in M with regular boundary dQ and a ET(n) is any local 
section (defined in an open subset containing Q). One defines then the "first variation" 
5S of S by considering homotopic variations rj = 5a G X U (' K ) with fixed values at the 
boundary dQ, x u i n ) being the space of vertical vectorfields of it. The critical sections 
are those sections along which 8S vanishes for any rj G X u i n ) with fixed values at the 
boundary. They are characterised by the equation: 

(A)*[e(£ o )] = 0, (1.3) 

which is the Euler-Lagrange equation. Here [e(£ )] is Euler-Lagrange morphism, a global 
bundle morphism e(£ ) : J 2 ^ ~^ A m (T*M) <8> V*n, where V*n is the dual bundle of the 
vertical bundle Vn, locally defined by: 

e(C ) = e A (C )ds ® d<t> A . (1.4) 

The first variation of Co is in fact globally defined through a further global bundle 
morphism f{C ) : J l n -> A m_1 (T*M) <8> V^*7r, locally expressed by: 

(j V) */(£„) = (p^K^o) o J V) ds, ® , (1.5) 

where ds^ is the surface element of dQ, defined so that ds, A dx^ = (— l) M <is. 
The following holds for T£ ; 

(jV)*[T(£o)] = (iV)* (e(C )\ V ) + (jV)*d(/(£ )k) , (1.6) 

for any local section a and any vertical vectorfield rj which projects onto a. Here and in 
the sequel (|) denotes standard dualitiy between forms and vectorfields. Equation (1.6) 
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is the global first variation formula of Co- The integral counterpart of equation (1.6) is 
thence the following 



6S= [ (e(C )\v)+ I (/0Co)h> 
Jn Jan 



(1.7) 



from which the Euler-Lagrange equations (1.3) follow after imposing the appropriate 
boundary condition r)\gn = 0. 

In order to study the stability properties of critical sections, i.e. of the solutions of 
the Euler-Lagrange equation (1.3), one considers next the variation of the action under 
second-order deformations of o. The second variation of S is then given by (see, e.g., 
[5]): 



S 2 S = - 



(e(£o)\p)+[ Hess £o (j 1 ?7 )+ f (f(C )\p) 
2 Un Jn Jan 

where p = 5 2 o~ denotes the second variation of o and the n-form Hess£ (j 1 7r) is the 
Hessian of C . Locally, equation (1.8) reads as follows: 



5 2 S = 



2 
+ 



dC A dC A 

P)A,A P + AAA P f 



d 2 C | L > 



d(p A d(p B 



d<P A , 
tTV + 2 



ds 



d 2 c, 



v A v B ,+ 



d 2 C 



d(j> A ^ B v 



°—ri A ri B 



ds 



.(1.9) 



A second order equation for n along critical sections, the (generalised) Jacobi equation, 
can be obtained from (1.8) by suitable integrations by parts on Hess. In fact, as discussed 
e.g. in [5], equation (1.8) can be conveniently rewritten as follows 



S 2 S 



1 



n JU x ' J dn 



1.10) 



where f(C )(rj,p) is a new boundary term depending both on r) and p and the term 
J&C£ (j 2 r}) is locally given by 

Jac £o (fv) = Mfv) d(p A , 



with 



d f dC B dC B \ ( , 10 



dx» \d(P A ^d(t) B 1 d4> A i p4> B v 
Equating JaU 2 v) to zero, i.e. setting 

Jac £o (A) = 0, (1.12) 
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gives rise to the standard form of the (generalised) Jacobi equation. 

Recalling that there is a natural bundle isomorphism X : V^n) — > ^(Vtt), locally 
defined by: 

(x\ y\ y\- r,\ V \) ~ (x\ y\ rf; y\, rj\) , (1.13) 
the first variation 5C : V(J l n) -> A n (T*M) of C , locally defined by: 

^^W^W;"- (L14) 

defines in fact a new Lagragian density in the bundle Vn by the prescription: 

d = 5C o (J)- 1 . (1.15) 

The new Lagrangian C\ is the first-order deformed Lagrangian. It is not difficult to prove 
the following: 

Theorem 1. Let (B,M,7v) be any fibered manifold and Co be any first order lagrangian 
density on it. Then the system formed by the Euler-Lagrange equation (1.3) and the 
Jacobi equation (1.12) of Cq is equivalent to the Euler-Lagrange equations of C\. 

This result, already used in [11] for the purpose of application to Riemannian Geometry, 
will not be deduced here. It will in fact follow from the analogous statement for second- 
order Lagrangians which will be discussed in the seubsequant Sections of this paper. 



2 Field Theory 

We discuss now the generalisation of these results, including the Jacobi equation as devel- 
oped in [5] , for general second-order field theories. For the sake of immediate understand- 
ing, the discussion will be given first in a coordinate language and using the analytical 
technique of truncated power expansions to express homotopic variations. Intrinsic geo- 
metric expressions will also be briefly discussed. The reader is referred to [5] for further 
details about notation. 



2.1 Series Development of the Action and Euler-Lagrange 
Equations 

Let us start by considering an action 

S = f C (x, 0, d(j), d 2 (p)ds, (2.1) 
Jn 

and let us consider an infinitesimal expansion of the action S around a section <p A . The 
infinitesimal deformations of <p A together with its derivatives <p A and (p A n U are given by 
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ev A + ^e* P A + o(e 3 ), 



4>\ = ^ A o) fi + ev\ + le 2 p\ + o(e 3 ), 
V = (0^) + o^V + o^ 3 ), 



(2.2) 



' w 2 

where e is a smallness parameter, while 77 = t] A 8a and p = p A &A are vertical vectorfields 
which correspond to the "classical" first and second variation 5(f) and 5 2 4> respectively. 
The variation of the action, at the second-order in e, is given by 



where 



and 



5S 
6 2 S 



= I 
Jn 



S{e) =S + e5S + -e 2 {5 2 S) + 



dC A 



dC 



5C 



d<j>\ 



A 



dC 



5<p 



„A 



d(f) A 
SC 



A 



\ P fi' fi^A P At" 



ds , 

ds + £2 ds , 
Jn 



(2.3) 



(2.4) 



Co = 



d 2 C 



d(fi A d(f) 
d 2 C 



H 2 r f> 2 r 
n A n B \ 2 n A n B +2 — - n A n B 



d(p A d(p B 



5- r/% + 2 



00V<^ 



2 £n 



2 £ o 



(2.5) 



is the Hessian of £0, which by later convenience we have denoted by £2- 

With the expansion above, we can now deal with the prob;em of extremality and 
stability of the action. Hereafter, the original Lagrangian £ will be simply denoted by £ 
to simplify notation. 

Let us then assume that (p A = <p A (x) is a section along which S has an extremum, 
i.e. the first variation of the action vanishes for all first-order deformations i] A . To this 
purpose it is convenient to perform an integration by parts and rewrite 5S as follows: 



SS= [ ^AV A ds+i 
Jn od) A Jan 



' 5C 



5C 



ds 



H J 



(2.6) 



where 
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6^ dx» \d(j) A J + dx^dx" \d^~) ~ dx^ \ 5 4> A J ' 

SC dC d ( dC \ 



5ct) A ^ d4> A ^ dx v \d(f> 



A 



are the Euler-Lagrange derivatives. 

As in the first -order case (see [11]), the local coordinate expression (2.6) corresponds 
to the intrinsic decomposition 

5S= I (e(C)\ V )+ f (/(£)|A), (2.8) 
Jn Jan x ' 

where e(C) : J 4 7r — > A m (T*M) <g> V*n is the Euler-Lagrange morphism and f(C) : J 3 n — > 
A m_1 (T*M) ® V*tx x is the canonical Poincare-Cartan morphism (uniquely existing for 
second-order theories, as is well known from the general theory; see [12], [13] and refer- 
ences quoted therein). 

The extremality of the action is classically described by the Euler-Lagrange equations 

^ = 0, (2.9) 



Scf) 

under the boundary conditions 



on 



an 



= 0, 

= o, 



(2.10) 



which just mean that the deformation of the local section O is fixed at the boundary, 
together with its first derivatives. The above conditions just concern the extremality of 
the action. They do not say anything about the stability of the solutions and to this 
purpose we must go to the next order of deformations. 



2.2 The Jacobi Equation 

A criterion for the stability of the Euler-Lagrange equations is obtained by looking at the 
sign of the second variation 5 2 S. If 5 2 S > we have a minimum; if S 2 S < we have a 
maximum; if 5 2 S = we might have a degenerate critical point. Assuming that we are 
"on shell", i.e. that the Euler-Lagrange equations hold along the section O , then S 2 S 
reduces to 

In order for 5 2 S to have a definite sign we impose the further boundary conditions 
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p 



an 



an 



0, 
0, 



(2.12) 



i.e., we require that 4>o is strongly fixed at the boundary up to order two. Under these 
stronger conditions 5 2 S reduces to 



S 2 S 



= [ 

Jo. 



Li ds . 



(2.13) 



As already discussed in [5] the Jacobi equation can now be obtained by a suitable 
integration by parts. Let us however remark that the quadratic form appearing in (2.13) 
can be changed by adding surface terms which, due to conditions (2.10) vanish at the 
integration boundary. Therefore, the form of the Jacobi equation to which one arrives 
strongly depends on the terms added at the boundary. One criteria would be to add 
surface terms which do not change the dependence on (77, dr], d 2 rj). For this let us rewrite 
£ 2 in the form 



d 2 C A B d 2 C A B d 2 C A B 

C2 = d(f> A d(f) B 71 71 + WW^ V 71 " + <9<W V 71 71 

d jCs a b d jCs a b OH a b 

+ ww v ^ + WW» v " v v + WW^x 71 1x71 uX 

d 2 C d 2 C d 2 C 

(2.14) 

The only terms which can be integrated by parts are the fourth, fifth, seventh and 
eighth terms. The four possible outcomes are summarised hereafter: 



L-0 = 



d ( dC d ( 8C 



d(p A \d(p B dx^ \d(p B . 



rj A rf + 



d 2 C 



d 2 C 



KB ^ A A 



v A v\ 



p,J 



+2 



d 2 c 



2 jO a b o d 2 jC, a b 

TA~^lB~ ^ ^ u + 1 + p., A ajB ' ^ M ^ ^ 



pv 



d 2 C 



pu^V \p 



B V B \p + 



d 



d 2 c 



dx» \d(t> A d(j) B ^ 71 J ' 



d 2 C 



V A V B +[2 
d 2 C 



d(p A d(p B 



+ 2 



d 2 C 



d 



d 2 c 



d(p A d(j) B dx 
d 2 C 



v A v B p 



d(f) A dcf) B 



[IV 



A »d<f> B u/ , 

f> 2 r 

n A n B 12 n A n B 



d 2 c 



d<p A 



Tb—V puV x P + 



d 



d 2 c 



dx^ yd^^l 



d<P A »d<P B vX 

jrV A v B )l , 
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C 2 * - M^rfrf+U-** - — ( —^—\\ ri A r) B 

2 ~ d<j) A d<f) B 1 1 \ d<p A d(j) B ^ d(P A dx» \d(f) B J 1 1 M 

d 2 C A B < d 2 C d 2 C \ A B 

+ d<j> A d<j> B ^ v ^ + \d^d^ v d<p A d<p B J v ^ " 

d 2 C d 2 C 
d ( d 2 C A B x 



+ dx» \d4> A d(t> B ^ v 7 ' 

£ ^ C r> A Tt B I 2 — T7 A T7 B +2 — ?7 A T7 B 

( d 2 c d ( d 2 c \\ A B 

I ( ) -A n 2 C A 

+ Me^^^ (2 ' 15) 

where = means equal modulo a total divergence. 

Apart from first-order divergencies, we can also add second-order divergencies to 
the Hessian. The only two possibilities which do not introduce derivatives higher than 
fourth-order ones are the following 



d 2 ( d 2 C A B N 



dx^dx v \d(f) A u d(f> B , 



d2 ( & C ^ ,A r , A d ( d 2 C \ au 



rj r] +4 — — - 7771— F-^- V V 



dx^dx" \d4> A ^d4> B u ) dx^ \d<f) A M <90 B v 

d 2 C d 2 C 

\d<f> A d<f> B ^ ' ' ) 

d 2 ( d 2 C \ A R „ d ( d 2 C d 2 C 



dx^dx" \d<j) A d<f) B 11 dx» \d(p A d<j) B ^ d(p B d(p A ^ . 
( d 2 C d 2 C \ d 2 C 

+ \ww;„ + ww; v ) vAr]B ^ +2 ww;/^ Bv - (2,16) 

A simpler possibility is to combine all the above integration by parts in a single 
decomposition displaying better properties. The result similar to the analogous one for 
first-order theories (see [5]) is given by 
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d ( 5C 



d<t> A \5(p B J 

( 2 C 

+ 
+ 



d 2 C 



d 2 C 



d 2 c d 2 c 



—\ v A v B 

ARAB \ 'I n'l v 



-jv A v\ 



V A uV B uX + 



d 2 c 



did 



5C 



,A R 



d<f> A p„d<f> B Xp 



n A n B 
'/ pu '/ Xp 



d 2 c 



-A B 



dx» \d<p A \5<j) B fl/ 



d 2 c 



pV 



d<P A ,d<p B uX 



— V A uV B \ 



(2.17) 



However, other terms can be added to the Hessian changing the dependence on only 
(77, dr], d 2r q). An interesting one is the following further expression: 



d 2 C 



B 



d 2 C 



KAQAB 



8(j) A d4> B , 



d 



d 2 c 



Vx + 
d 2 C 



d 2 C 



d<f> A d<j) B Xp 



dx^ \d(P A ^d(j) Br] + <90^<90 B A 77 A + 



V x P 
d 2 C 



B 



d 2 



d 2 C 



dx^dx v 



rj + 



d 2 C 



^ A pM B x 



V B x + 



Xp 

d 2 C 



pv^Y Xp 



d 



dx^ 1 
d 



d 2 C 



d 2 C 



d 2 c 



d<t> A M B x 



d 2 C 



B 



d 2 C 



dx» yd^d^ 71 1 dct> A ^ B x 



V B x + 



d 2 C 



A> v A 



+ 



pV^T J pv 

d 2 C B d 2 C B 



pv^ 



pv u V X 



d<P A pM B xp 
92C v B ) v A 

'I Xp I 'I v 



d^ A pM B xp 



(2.18) 



In all the above cases the surface terms cancel at the integration boundary, therefore, 
different Jacobi equations would be obtained. The equation ensuing from (2.18), namely 



d 2 C 



d 2 c 



R 



d 2 C 



d(p A d(j) 



B 71 + 8(p A d(j) B x 11 X+ £>shAf)AB_ V Ap 



d 



dx^ \d<fi A p d(f> 
d 2 ( 8 2 C 



d 2 C B d 2 C u 



d(f> A d(f> B Xp 
d 2 C 



p u V X 



rj B + 



d 2 C 



dx»dx v \d(j) A ^d(j) B 1 d<p A uu d(P B 



V B x + 



d 2 c 



pv u V x 



d^p^xp 



V B x P 



0, (2.19) 



is the "standard Jacobi equation" for C. It is the equation (0.2) we already mentioned in 
the Introduction 
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3 Hierarchical Structure of the Deformed Lagrangian 



Now we introduce a hierarchical structure associated to the deformations of the La- 
grangian. 

3.1 The Deformed Lagrangian 

Let us again consider an infinitesimal deformation of the field <f> A and its derivatives 4> A „ 
and (f) A ' up to second-oder in e. Then, the variation of the Lagrangian, at second 
order-order in e, is given by 

C( V )=C + eC 1 (7 ] ) + ±e 2 [C2(v)+Mp)\ > (3-1) 

where L\ is given by 

~ / \ 9C A dC A dC A 

CM = w v w; v ^w;, 11 (3 - 2) 

while Li is given by equation (2.5). 

The introduction of the above definitions is not only a matter of notational convenience 
since Co, C\ and £ 2 satisfy a series of remarkable identities. In fact, it can be checked 
that the Euler-Lagrange derivatives of L\ and £ 2 are given by 

5r) A S(f> A ' 
5rj A 8(j) A 

Moreover, it can be checked (see also [4]) that the equation 

5C 2 



(3.3) 



v °' (3 - 4) 

corresponds to the Jacobi equation (2.19) in its standard form. 

3.2 Unified description of the Euler— Lagrange and Jacobi equa- 
tions 

Let us now remark that in virtue of (3.3) the Euler-Lagrange equation and the Jacobi 
equation can be rewritten together as 



8Co 

5<P A 

Sr] A 



6C1 
5r] A 
6£i 



0, 



= 0. 



(3.5) 
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Therefore, we can consider the above system of equations as coming from a single vari- 
ational principle based on the deformed Lagrangian C\ where, however, the dynamical 
variables have been "doubled" to 0's and 77's. 

We can then introduce momenta canonically conjugated to 0's and 7/'s by means of 
standard prescriptions 



7T A - = 



Pa' 



n A 



w, 

dCi 

dp 
ggi 
dv A ^ 

dr\ A 



dx v 



d ( dCi 



dx v 



drj A 



and a canonical energy-momentum tensor by setting: 

W"„ = ir A » <P A V + N A <* <f> A Xv + pa" V a u + n A ^ V A Xu -6Z&. 



(3.6) 



(3.7) 



The explicit expression of Ti^v for our case is nor relevant here. What is really im- 
portant is the fact that, if Co does no depend on "space-time variables" x x , then TC^ U is 
a conserved quantity. Therefore, the allowed deformations are selected by the first-order 
condition 



dx» 



= 



(3.8) 



which can thence replace the second-order condition given by the Jacobi equation (2.19). 
Let us also remark that the Hessian matrix associated to £1 is given by 



a 2 d 



d 2 d 



\ 



(3.9) 



so that the following holds true 



|det[W(A)]| = [det[W(C )]] 2 , (3.10) 

which ensures that the regularity of C\ depends only on the regularity of C . 

Suitable globalisations of this concept can be achieved by introducing a connection 
(see, e.g., [14], [15]). See also [16] and [17] for a discussion of "rf-invariance" of these 
energy-momentum tensors. 



14 



4 Examples and Applications 

Here we consider some simple Lagrangians which involve up to second-order derivatives 
in order to illustrate our general results. 

4.1 The Riewe Lagrangian 

The Riewe Lagrangian (see [18]) was introduced to describe a classically spinning particle. 
It is given by 



1 . 2 1 m 
-ma — - — 

2 H 2 u 

where q is a three-dimensional vector. The equtions of motion are given by 



_ = -,^ f + _^=0. (4.2) 
The solutions of the above equation of motion are given by 



q i = + v ^t + a 1 cos(u; t) + b l sin(w t) . (4.3) 

These solutions describe a helical motion along a ellipsoid of semiaxis a 1 and b % oriented 
in the direction of v l . 

The corresponding Jacobi equation is given by 

rf + i-^ = 0. (4.4) 

4.2 The Shad wick Lagrangian 

The Shadwick Lagrangian is given by 

£ = (0oo0n-0oi) , (4-5) 
where is a scalar field (see [19]). The field equations are given by 

5C 



vu =3 (0oo0n-0oi) =0- (4.6) 
The corresponding Jacobi equation is given by 

0ii »7oo + 0oo »7n - 2 0oi mi = ■ ( 4 - 7 ) 
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4.3 The Klein— Gordon equation 

The Klein-Gordon equation is obtained from the Lagrangian density 



A) = \ (<T h 4>u - m 2 2 ) , (4.8) 



2 

whose first-order and second-order deformed Lagrangians are given by 



£* = \ (<T^-™V) • (4-9) 



The Euler-Lagrange and Jacobi equations are given by 

6C1 



= (<ro in , ■ nro) 0, 

-(g^V»» + m 2 v) =0, (4.10) 



5<p 

respectively. Solutions to these equations are expressed by 

4>{X, t) = J e Wx+^+rnH) + e i(h-W()j dA . j 

r/(x, i) = I (fc+(fc) e ^x+VkWt) + e i(fcx-v^+m^)) ^ , (4.11) 

such that 

o (x, t) = 1 J Vk 2 + m 2 (f + (k)e l{kx+Vk2+mH) - f_{k)e l{ - kx -^ k ' 2+m2 ^) dk , 
^(x, t) = t j k t (f + (k)e l{kx+ " /W ^ t) + f4k)e t{kx ~ VWT ^ t) ) dk, 
T) (x,t) = i J Vk 2 + m 2 (h + {k) e W Q - h_(k) e ^x-^+ m H)^ dk ^ 
rn(x,t) = 1 J h {h + {k)e ii - kx+VW ^ t) + h^{k)e^ kx - Vk2+m ' 2t ^ dk . (4.12) 
The appropriate energy-momentum tensor is thence given by 

W„(A) = ^ ^ x 0a + Vu vx - % A , (4.13) 

and its time-time component 7i° (£i) is given by 

7A(£i) = <poVo + <fiiVi + m 2 <PV ■ (4.14) 
The associated energy (see also [6], [7]) is therefore given by 
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and an evaluation on-shell gives immediately 



(4.15) 



E(d) = 2(27r) 3 m 2 / (f + (k) h-(h) + f-(k) h+(k)) dk = constant 



(4.16) 



This condition is equivalent to the Jacobi equation. 
4.4 Geodesies in a Riemannian manifold 

A further example of application to standard Jacobi fields along geodesies of a Riemannian 
manifold (M,g) has been discussed in [11]. Here we can shortly summarise the results. 

Let Q be a n-dimensional manifold and (TQ, Q, tq) its tangent bundle. Let g be a 
Riemannian metric on Q. The geodesies of (Q,g) are those cuves 7 : R — > Q whose 
tangent vector 7 is parallel along 7, i.e., it satisfies V^,7 = 0; in local components 



The Jacobi fields of (M,g) are those vectorfields n = rfdi, defined along geodesies 7 
by the differential equation: 



where V 2 denotes the second-order covariant derivative along the curve 7 and Riem is 
the tri-linear mapping defining the Riemannian curvature of g. Jacobi fields generically 
define infinitesimal deformations of geodesies into families of nearby geodesies. According 
to [20], the metric g can be lifted to a metric g c on the manifold TQ, the "complete 
lift", as follows. Let g = gijdq l dq> in a local chart (U;q l ); then the corresponding local 
expression of g c in (TU; q l ,u l ) is 



0. 



(4.17) 



V 2 rj + Riem(?7, 7, 7) = , 



(4.18) 



g c = 2g ij 5u i dq j , 



(4.19) 



where 5u l denotes the following 



5u l = du { + T^fc u m dq k . 
For any function / : Q — > R a new function df : TQ — > R is defined by setting 



(4.20) 




(4.21) 



g c = (dgij) du l dv? + 2 g^ du 1 dq 



(4.22) 



i.e., the (2n x In) matrix of g 



is 
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f-C* f). (4.23) 



The following holds true: 



Theorem 2. Let (Q,g) be a Riemannian manifold. Then the system formed by the 
geodesic equation of g in Q and the Jacobi equation associated to g in TQ is the geodesic 
equation in TQ of the complete lift metric g c . Therefore this system follows from a 
variational principle on TQ based on the energy functional defined by the lifted metric g c . 

We see immediately that this theorem is nothing but a simple consequence of our general 
results. In fact, the energy functional of g is based on the Lagrangian: 

C = \g l3 u l u\ (4.24) 
whose associated first-order deformation Lagrangian is thence given by 

A = \ {dkQij) u l v? V k + 9ij u* if . (4.25) 
Using V(<?) = this becomes immediately 

U=9a [v l + r mk u m V k ] v? . (4.26) 
Then C± is in fact the energy Lagrangian of the lifted metric g c = 2 gijSrfdqi . 



5 Conclusions 

We have presented here a natural and direct generalisation of the Jacobi equation to the 
case of second-order Lagrangians, which are important not only for the sake of complete- 
ness but also and specially in view of applications to relativistic field theories (whereby 
gravitational Lagrangians depend effectively of second-order derivatives of a metric). 
Concrete applications to relativistic field theories will form the subject of forthcoming 
investigations. 



Acknowledgements 

One of us (B. C.) acknowledges support of GNSAGA-CNR and of MURST (Nat. Res. 
Proj. "Geometria delle Varieta Differenziabili"). 

One os us (M. F.) acknowledges support of GNFM-CNR and of MURST (Nat. Res. 
Proj. "Metodi Geometrici e Probabilistici in Fisica Matematica"). 

One of us (V. T.) has been partially supported by a Visiting Professorship of the 
Italian GNFM-CNR and of the International Centre for Theoretical Physics, Trieste. 



18 



References 



[1] C. Lanczos, The Variational Principles of Mechanics, 4th Edition (University of 
Toronto Press, Toronto, 1970). 

[2] W. Klingenberger, Riemannian Geometry (W. de Gruiter, Berlin, 1982). 

[3] H. Rund, The Hamilton- J acobi Theory in the Calculus of Variations (Van Nostrand, 
Princeton, 1966). 

[4] H. Taub, Stability of general relativistic gaseous masses and variational principles, 
Commun. Math. Phys. 15, 235 (1969). 

[5] B. Casciaro and M. Francaviglia, Covariant second variation of first order La- 
grangians on fibered manifold. I. Generalized Jacobi fields, Rend. Mat. Appl. 16, 
233 (1996). 

[6] O. Amici, B. Casciaro and M. Francaviglia, Covariant second variation for first order 
Lagrangians on fibered manifolds. II. Generalized curvature and Bianchi identities, 
Rend. Mat. Appl. 16, 637 (1996). 

[7] O. Amici, B. Casciaro and M. Francaviglia, Second variation and generalized Jacobi 
equations for curvature invariants, Atti Accad. Peloritana dei Pericolanto 74, 73 
(1996). 

[8] O. Amici, B. Casciaro and M. Francaviglia, The second variation for non-linear 
gravitational Lagrangians, Atti Accad. Sci. Torino CI. Sci. Fis. Mat. Natur. 130, 
263 (1996). 

[9] P. Foulon, Geometrie des equations differentielles du second ordre, Ann. Inst. H. 
Poincare 45, 1 (1986). 

[10] E. Martinez, J. F. Carinena and W. Sarlet, Derivations of differential forms along 
the tangent bundle projection. II, Diff. Geom. Appl. 3, 1 (1993). 

[11] B. Casciaro and M. Francaviglia, A new variational characterization of Jacobi fields 
along geodesies, Ann. Mat. Pura Appl. 172, 219 (1997). 

[12] M. Ferraris, Fibered connections and global Poincare-Cartan forms in higher-order 
calculus of variations, in Proceedings of the Conference on Differential Geometry 
and applications, Nove Mesto na Morave, 1983, ed. D. Krupka (J. E. Putkyne 
University, Brno, 1984), pp. 61-91. 

[13] I. Kolar, A geometrical version of the higher order Hamiltonian formalism in fibered 
manifolds, J. Geom. Phys. 1, 127 (1984). 

[14] R. Ragioneri and R. Ricci, The Legendre transformation, Boll. Un. Mat. Italiana B 
18, 119 (1977). 



19 



[15] M. Ferraris and M. Francaviglia, On the globalization of Lagrangian and Hamilto- 
nian formalisms in higher order mechanics, in Proceedings of the IUTAM-ISIMM 
Symposium in Modern Developments in Analytical Mechanics, Torino, 7-11 June 
1982, eds. S. Benenti, M. Francaviglia and A. Lichnerowicz (Tecnoprint, Bologna, 
1983), pp. 109-125. 

[16] V. Tapia, N. Ferraris and M. Francaviglia, Identically vanishing field equations and 
d-invariance in field theory, Nuovo Cimento B 103, 435 (1989). 

[17] M. Ferraris, M. Francaviglia and V. Tapia, Global d-invariance in field theory, J. 
Phys. A: Math. Gen. 26, 433 (1993). 

[18] F. Riewe, Generalized mechanics of a spinning particle, Lett. Nuovo Cimento 1, 807 
(1971). 

[19] R. Shadwick, The Hamiltonian formulation of regular rth-order Lagrangian field 
theories, Lett. Math. Phys. 6, 409 (1982). 

[20] K. Yano and S. Ishihara, Tangent and Cotangent Bundles (M. Dekker, New York, 
1973). 



20 



